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ABSTRACT 

Context. The supernova remnant (SNR) RX J1713. 7-3946 (also known as G347. 3-0.5) is part of the class of remnants dominated by 
synchrotron emission in X-rays. It is also one of the few shell-type SNRs observed at TeV energies allowing to investigate particle 
acceleration at SNRs shock. 

Aims. Our goal is to compare spatial and spectral properties of the remnant in X- and y-rays to understand the nature of the TeV 
emission. This requires to study the remnant at the same spatial scale at both energies. To complement the non-thermal spectrum of 
the remnant, we attempt to provide a reliable estimate for the radio flux density. 

Methods. In radio, we revisited ATCA data and used HI and mid-infrared observations to disentangle the thermal from the non- 
thermal emission. In X-rays, we produced a new mosaic of the remnant and degraded the spatial resolution of the X-ray data to 
the resolution of the HESS instrument to perform spatially resolved spectroscopy at the same spatial scale in X- and y-rays. Radial 
profiles were obtained to investigate the extension of the emission at both energies. 

Results. We found that part of the radio emission within the SNR contours is thermal in nature. Taking this into account, we provide 
new lower and upper limits for the integrated synchrotron flux of the remnant at 1.4 GHz of 22 ly and 26 Jy respectively. In X-rays, 
we obtained the first full coverage of RX J1713. 7-3946 with XMM-Newton. The spatial variation of the photon index seen at small 
scale in X-rays is smeared out at HESS resolution. A non-linear correlation between the X- and y-ray fluxes of the type Fx °c F^' 41 is 
found. If the flux variation are mainly due to density variation around the remnant then a leptonic model can more easily reproduce 
the observed X/y-ray correlation. In some angular sectors, radial profiles indicate that the bulk of the X-ray emission comes more 
from the inside of the remnant than in y-rays. 

Key words. ISM: supernova remnants - Supernovae : individuals : RX J1713.7-3946 - Acceleration of particles 



[ 1. Introduction RX J1713.7-3946 was first dis covered in X-rays with the 

. ROSAT all-sky survey in 1996 dPfeffermann & Aschenbachl 

*_h ' Il996l) . The remnant is close to the galactic plane and its distance 

Supernova remnants (SNRs) have long been believed to be ac- is controversial. Using ASCA observations and the measurement 

celerators of cosmic rays at least up to the knee (~ 10 15 eV). of the column density to ward the source , it was first proposed to 

Evidence that electrons are indeed accelerated in SNRs is found be ata distance of 1 kpc (|Koyama et alj[1997|). Then|Slane etap 

both in radio and X-rays through synchrotron emission. The de- (1922) derived a distance of 6 kpc based on a possible associ- 

tection of TeV emission from a SNR is an evidence of proton ation of the remnant with molecular clouds and the HII region 

acceleration if the y-rays come from the interaction of acceler- G347.6 1 1 +0.204. The comparison of the X-ray and HI absorb- 

ated protons with the ambient matter (hadronic model). However in § column densities dCassam-Chenai et al|2004|) as well as CO 

the y-rays can also be produced via Inverse Compton scattering observations of a molecular cloud interacting with the remnant 

(IC) of accelerated TeV electrons off ambient photons (leptonic dFukui et alj2QQ3|) both suggest a closer distance of 1 kpc (value 

model) that could be either infrared surrounding emission or cos- adopted in this paper). In this case, the remnant is about 20 pc 

mic microwave background. The SNRs interacting with dense in diameter (70 arcmin on the sky) and could be associated with 

molecular clouds are good candidates to detect y-ray emission the supernova that exploded in AD 393 in the tail of the constel- 

from the hadronic mechanism as the high density of the cloud lation Scorpius dWan S et alj|l997|). The remnant would then be 

provides a large amount of targets for the accelerated protons, about 1600 years old; 

TheSNRRXJ1713.7-3946(alsoknownasG347.3-0.5)isoneof In me radio band|Slane et al.| 41999b showed for the first time 

those candidates as it is interacting with a dense molecular cloud an ima g e of RX J 1 7 1 3 .7-3946, bas e d on MOS T observations 

in the northwest (NW) and in the southwest (SW) of the remnant f 843 MHz. Later on, |Ellison etal.| 4200JJ) and |Lazendic et afl 

dFukui et alJl200l ICassam-Chenai et alJl20M iMoriguchi et all (with a recalibrated and improved image) reported ATCA 



|2QQ5|) . radio observations of this SNR at 1.4 GHz (left-hand panel of 
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Fig. [5]). The SNR appears as a faint nebula (~ 1° in diame- 
ter) with many short, curved features, the brightest of which 
are the two bright arcs visible on the northwest corner. These 
arcs coincide with the edges of the brightest X-ray emission. 
Interestingly, those arcs are observed to be in the vicinity of 
the HII region G347.6 1 1 +0.204 which is located at 6.6 kpc 
dRusseill 12001 IConti & Crowtherl I2004T) . An almost complete 
ring of weak emission, about 30 arcmin in size, can be also de- 
tected near the center of this extended SNR. An accurate esti- 
mate of the integrated radio flux density of RX J1713. 7-3946 is 
a pending problem. This is due to the intrinsic faintness of this 
SNR, the possible mixing of the SNR synchrotron emission with 
the thermal emission of the nearby HII region as well as the lim- 
ited quality of the available data. 

In X-rays the emission is dominated by a non-thermal 
continuum and no emission lines have b een observed so far 
dKovama et al l fl997| ISlaneetalJ [19991: iPannuti et all l200l 
Cassam-Chen ai et al.ll2004 ). This non detection can set an upper 



limit to the ionization age. Assuming tsNR = 1600 yrs an upper 
limit on the d ensity of the ambient medi um of 0.02 crrr 3 is 
then derived dCassam-Chenai' et alj l2004t) . The study of the 
re mnant with XMM-Newton at sm all spatial scale car ried out 
bv lCassam-Chenai et al.1 (120041) and lHiraga etal.1 (120051) clearly 
showed a spatial variation of the photon in dex. It varies from 1 .9 
to 2.6 with a mean statistical error of 4% dCassam-Chenai' et al.l 
2004). A correlation between the X-ray flux and the photon 
index is also observed. Those results were obtained using an 
adaptive grid to have approximately the same number of counts 
in each pixel grid region where a spectrum is then extracted. On 
the bright regions (the north-west rim), the size of the grid pixel 
is typically 0.05°. 

In y-rays, RX J1713. 7-3946 has been detec te d by 
the CANGAROO collaboration dMuraishi et all l2000t 
Enomoto et al.l 120021) . It was then definitely established by 
the HESS telescopes which provided the first spa tially resolved 
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y-ray image of the remnant dAharonian et al . 2004). The overall 
morphology is very similar to that in X-rays (same western 
shell) and the brightest spots (located to the north-west) are 
coincident. As the remnant is extended and is the brightest 
shell-type SNR seen in y-rays, its emission can be studied in 
detail. With the good spatial resolution of the HESS telescopes 
(~ 0.1 °) it has been possibl e to carry out a spatial ly resolved 
spectral study at large scale ( Aharoni an et al.l 1 20 06, hereafter 
AH06). No evidence for spatial variations of the photon index 
from region to region (mean value of 2.09 and a standard 
deviation of 0.07) was found down to the precision of the 
measurements (mean l<x statistical error of ± 0.08). Moreover 
there does not seem to be any correlation between the y-ray flux 
and the photon index (Fig. 14 of AH06). 

Those y-ray results are very different from what has been found 
at small scale in X-rays. However it is important to note that 
the X- and y-rays results were not obtained looking at the same 
spatial scale (the spectra were not extracted using the same size 
of extraction regions). Moreover the Point Spread Function of 
both instruments is very different. Those differences hamper a 
comparison of the results at both energies. The major issue here 
is whether the spectral parameters are really different in X- and 
y-rays or if this difference is introduced by the fact that both 
studies are looking at different spatial scales. 
To investigate this question we carried out a detailed compar- 
ison of the remnant in X- and y-rays using the same extraction 
regions and taking into account the different spatial resolution of 
the two instruments. We compared both the spectral properties 
of the remnant and the morphology using radial profiles in X- 



and y-rays. 

Concerning the radio flux estimate, we investigated which 
part of the radio emission is synchrotron emission in origin 
(related to the remnant) or thermal emission (possibly related to 
the nearby HII region). We then calculated the integrated flux 
density re-analysing Lazendic et al. (2004)'s database. 
In Sect. [2] we present the different sets of data used in this study 
including the new XMM-Newton observations that complete the 
mosaic of the remnant. Section[3]presents the methods used for 
the processing of the X-ray data. In Sect. |4] the estimate of the 
integrated radio flux of the remnant and the results of the X- 
and y-ray spectral and morphological comparison are shown. 
In Sect. [5] we discuss this comparison in the framework of an 
hadronic and leptonic model. The implications of the new radio 
flux for the multi-wavelength emission models of the remnant 
are also discussed. 



2. Data 

2.1. Radio data 

In this w ork we re-analyzed ATCA observations carried out at 
1.4 GHz dLazendic et al.ll2004l) following different paths to esti- 
mate the radio flux density of RX J1713. 7-3946. The details to 
carry out an accurate calculation are discussed in Sect. 14. II 



2.2. Infrared data 

Infrared observations in and around the remnant can help dis- 
entangling the thermal from the non thermal radio emission 
which is crucial to estimate the radio flux of the whole rem- 
nant. We used a method ba sed on a color-color criteria pro- 
posed by iReach et al.1 d2006) using Spitzer data at 3.6, 4.5, 5.8, 
and 8 fim from the Galactic Legacy Infrared Mid -Plane Survey 
Extraordinaire (GLIMPSE, IBeniamin etaT]|2003l) . The spatial 
resolution of the survey is about 2 arcsec. More details can be 
found in Sect. ED 



2.3. XMM-Newton data 

The first set of observations of the SN R RX J 17 13.7-3946 (pre- 
sented in Cass am-Chenai et a l. 2004) covered almost the entire 
remnant. New observations of the south, east and west regions 
completed the mosaic of the remnant. With a total of 1 1 point- 
ings (see Table [D we can now provide the first full mapping of 
RX J1713. 7-3946 with XMM-Newton (see Fig. [T}. 
To clean proton flare contamination in the event files, we built 
a histogram of counts/seconds in the 10-12 keV band for the 
MOS cameras (12-14 keV for PN). Then we fit a gaussian dis- 
tribution upon the histogram and retain, in the observation, only 
the time intervals whe re the count rates are within a 3cr range 
dPratt & Arnaudll2002l) . Table [T] provides the list of the remain- 
ing exposure times after flare screening. For the PN instrument 
only 6 pointings were used. 



2.4. HESS data 

In the HESS study (AH06) spectra were extracted from 14 
square regions (0.26° length, see Fig. |2]l covering the whole 
remnant. The data used for the spectral analysis have a cut on 
the minimum size of image of 80 photo-electrons resulting on a 
mean spatial resolution of 0.12° (68% containment radius). For 
the morphological analysis (radial profile) the cut used in HESS 
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data is at 200 photo-electrons resulting on a better spatial res- 
olution of 0.08°. The detection efficiency varies over the 5° of 
the full field of view but it is almost constant at the scale of the 
remnant which is ~ 1° wide (it varies by only 5% between the 
centre and the edge of the object). All the data used in our study 
are taken from AH06 and we did not reprocess any TeV data. 
For the comparison of the spectral properties we used the y-ray 
results presented in Table 2 of AH06 and for the comparison of 
radial profile we used the data presented in Fig. 16. 



3. X-ray processing 

3.1. Mosaic construction 

The mosaic is built in counts and an adaptive smoothing is ap- 
plied such that the signal-to-noise ratio is at least 10. The in- 
strumental backgro und is derived from a compilation of blank 
sky observations (Carter & Read 120071) . renormalized in the 
10-12 keV energy band for the MOS cameras (12-14 keV 
for PN) and subtracted from each image. To have units in 
photons/cm 2 /s/pixel instead of counts/s/pixel, the exposure map 
of each observation is multiplied by the average effective area 
in the energy band (assuming the same spectrum over the field 
of view). Then a mosaic of those exposure maps is built and 
smoothed in the same way. The final image is the division of the 
counts mosaic by the exposure map mosaic. 
The resulting image which is the sum of the MOS and available 
PN data after flare screening is presented in Fig.Q] The morphol- 
ogy of RX J17 13.7-3946 can be decomposed in two main kinds 
of structures : diffuse emission present over all the remnant and 
bright filaments particularly visible in the west and north of the 
remnant (see Fig. [TJ. Thanks to the high sensitivity of XMM- 
Newton, we can clearly see the faint emission in the recent ob- 
servations of the regions south, east and north. In particular in 
the northern region of the remnant we distinctly see a straight 
edge that is not an artifact due to a CCD gap or any instrumental 
effect. Also it is not due to an X-ray absorption along the line 
of sight as this straight edge remains visible on the 4.5-7.5 keV 
image (in this energy band the absorption is weak). Simply the 
emission seems fainter there. Above that edge, we see a structure 
(in blue-green) that seems to be the continuity of the shock. 

The estimate of the astrophysical background is not sim- 
ple in RX J17 13.7-3946 as it seems to vary around the rem- 
nant. However in order to have a rough approximation of the 
background level we extracted the flux outside a circle of 0.56° 
radius centered on the remnant (arj2ooo = 17hl3m46s, <5j2ooo = 
-39°44'56"). We then subtracted the mean value of this flux 
(4.5 x 10~ 6 photons/cm 2 /s/arcmin 2 ) for the morphological study. 
The small structures of the SNR are not affected by this back- 
ground subtraction. 



3.2. Spectral extraction method 

With its good spatial resolution, the XMM-Newton telescope 
can carry out spectral study at small scale whereas in y-rays 
the spectral analysis is done at larger scale due to the com- 
paratively lower spatial resolution of the HESS telescopes. To 
address this problem, we took into account the different Point 
Spread Functions and the variation of the detection efficiency 
across the field of view of the two instruments. 

In the case of XMM-Newton, the detection efficiency of the 
MOS and PN cameras can drop 35% from the centre to the edge 
of our 0.26° extraction regions. For the HESS telescopes the de- 
tection efficiency is almost constant to this size (see Sect 12.41 ). 



Table 1. XMM-Newton observations used in this paper. The total 
and good columns represent the exposure time before and after 
flare screening. 







Exposure (ks) 






MOS 


PN 


Obsld 


Observation Date 


Total 


Good 


Good 


0093670101 (NE) 


2001 September 5 


15.3 


1.8 





0093670201 (NW) 


2001 September 5 


15.3 


6.7 





0093670301 (SW) 


2001 September 8 


15.3 


15.2 


10.0 


0093670401 (SE) 


2002 March 14 


14.1 


11.6 


5.1 


0093670501 (CE) 


2001 March 2 


13.8 


13.0 


6.5 


0207300201 (CE) 


2004 February 22 


31.5 


12.4 


0. 


0203470401 (NE) 


2004 March 25 


17.0 


16.1 


6.7 


0203470501 (NW) 


2004 March 25 


18.0 


13.1 


9.7 


0502080101 (E) 


2007 September 15 


34.6 


5.8 





0502080301 (W) 


2007 October 3 


8.9 


2.8 





0551030101 (S) 


2008 September 27 


24.9 


24.5 


20.8 




12:00.0 17:10:00.0 



Right ascension 



Fig. 1. EPIC MOS plus PN image in the 0.5-4.5 keV band. The 
units are ph/cm 2 /s/arcmin 2 and the scale is square root. The im- 
age was adaptively smoothed to a signal-to-noise ratio of 10. 
The four ellipses show the regions used to estimate the local as- 
trophysical background for the spectral analysis. 



When extracting an X-ray spectrum in this large region we want 
all the events to contribute with the same weight to the spectrum. 
T o address this proble m, we used the weight method described 
in lArnaud et alJ (2001) where each event is corrected for its ef- 
ficiency loss as a function of its position on the camera and its 
energy. 

We also have to take into account the different size of the 
Point Spread Function of both instruments. For the spectral study 
of this SNR, the mean spatial resolution of the HESS instru- 
ment is 0.12° (68% containment radius) which is comparable to 
the size of the extraction region (0.26°) whereas with the XMM- 
Newton observatory, the spatial resolution is about 7 arcsecs. 

We assumed that the Point Spread Function of the HESS 
telescopes is a gaussian of <x = 0.0795° (corresponding to the 
HESS 68% containment radius of 0.12°) and that in comparison 
the Point Spread Function of XMM-Newton is negligible. To de- 
crease the spatial resolution of the X-ray data, we redistributed 
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Fig. 2. Spatial X-ray contribution to each HESS region (as defined in lAharonian et all I (120061) 1 We see that due to the size of the 
PSF comparable with the size of the extraction regions, many events outside of the region contribute to the spectrum. To save space, 
each frame shows the contributions to three or four separate regions at once. The linear scale is in ph/cm 2 /s/arcmin 2 . The regions 1 
to 14 defined here are used later for the spectral analysis in Table[2] Fig. [7] [8] and [9] 



randomly the position of each X-ray event according to the gaus- 
sian probability density function. 

For the purpose of spectral comparison with the y- 
rays, we removed the two following very bright point-like 
sources in the X-ray data : 1WGA J 17 13.4-3949 which 
is argued to be t he Central Compact Object of the SNR 
(Slan eetal.1 119991: lLazendic et all 120031: ICassam-Chenai et al.1 
12004 and 1WGA J1714.4-3945 w hich is associated with a star 
dPfeffermann & AschenbachfT9 96). 



3.3. X-ray background spectra 

The instrumental b ackground spectrum is derived from blank 
sky observations dCarter & Read! [2007b in the same detector 
area, renormalized in the 10-12 keV energy band for the MOS 
cameras (12-14 keV for PN) and processed with the same 
method as the observations (see Sect. 13.21 ). 

The subtraction of the local astrophysical background, is 
not a simple problem in RX J1713.7-3946. On the one hand, 
as the remnant is ~ 1° in size and close to the galactic plane, 
the local astrophysical emission can vary around the remnant. 
On the other hand, there are only few pointings that allow us 
to estimate this background. We extracted background spectra 
outside of the SNR in 4 regions where the statistics was suffi- 
cient (see the 4 ellipses in Fig. Q3. In this astrophysical back- 
ground, no emission lines are seen but the statistics is low. For 
comparing the spectral properties of the background between 
the different regions, we modeled the spectra using a simple 
absorbed power law model. The parameters of this model are 
similar in the northern, north-western and south-western regions 
i.e. : A^h — 0.4 x 10 22 crrT 2 , index ^1.9 and a normalization 
norm 1.4 x 10~ 3 keV _1 cm _2 s _1 at 1 keV (renormalized to 
the area of the large HESS extraction regions). Note that in 
the south-eastern region, the normalization parameter is 40 % 
smaller. As all the extracted background spectra have low statis- 
tics, we decided to use a power law model with the parameters 
listed above for the astrophysical background instead of a spec- 
trum. We used the same A^h, index and normalization parameters 
that of the northern, north-western and south-western regions for 
all the extraction regions. The impact of such a choice for the 
south-eastern region, where the normalization is smaller, is dis- 
cussed later in Sect. 13.41 



3.4. Spectral modeling 

The X-ray emission of the remnant is dominated by a non- 
ther mal continuum and no emission lines have been detected so 
far dKovama et al.1 1 19971: ISlane et all 1 19991: iPannuti et al J 120031: 
ICassam-Chenai' et al.ll2004l) . In our study, all spectra are well 
described by an absorbed power law model as is illustrated in 
Fig. [3] In some observations, the PN instrument exposure time 
was null after flare screening. To have an homogeneous cover- 
age of the remnant we kept only the MOS data. Whenever it was 
possible we fitted independently, as a test, the spectrum extracted 
from the MOS 1 &2 and PN instruments and compared the result- 
ing best-fit parameters. They agreed within the statistical errors 
bars. 

All the data was fitted using unbinned spectra with the C-statistic 
implemented in Xspec (vl2.3.1). Binning was used for graphi- 
cal purposes only and fixed at 3cr for all spectra. We have fit- 
ted all the data from 0.8 keV to 10 keV. The best-fit parame- 
ters for the HESS regions are listed in Table [2] As discussed in 
Sect. 13. 31 the local astrophysical background in the south-eastern 
region seems weaker than the average background used for all 
the regions. Therefore, we studied the impact of those different 
backgrounds for the faintest region (the most dependant to back- 
ground subtraction) which is located in the SE : region 12. When 
using a value of the normalization parameter 40 % smaller for 
the background model, the new best-fit parameters do not vary 
much. While the hydrogen column and the index change within 
the errors bars, the integrated X-ray flux varies by less than 10%. 
As the variations are small and do not impact the results of our 
study, we kept the same astrophysical background for all the re- 
gions. 



3.5. Radial profile method 

As we now have a full coverage of the remnant, we can properly 
compare the radial profiles of the emission of the remnant in X- 
and y-rays. To match the HESS spatial resolution we smoothed 
the X-ray image with a gaussian of cr=0.053° (corresponding 
to r6g%=0.08°). We then extracted a flux per unit solid angle as 
function of the distance to the center on the smoothed X-ray im- 
ages (1-2 keV and 2-4.5 keV band) in eight sectors as shown in 
Fig. [4] To compare the radial profiles at both wavelengths, the 
X-ray profiles were scaled by a unique normalisation factor cal- 
culated as the ratio of the total number of counts in y-rays over 
the total flux in X-rays on the whole remnant. In the last two 
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10 FT 




1 2 5 

Energy (keV) 

Fig. 3. Best-fit X-ray spectrum from region 3 (Top panel) and 
region 12 (Bottom panel) with an absorbed power law model 
(MOS1 spectrum is in black and MOS2 in red). The region 3, 
located on the north-west, is the brightest region of the rem- 
nant. It has a steep spectrum (2.35) and has a high absorption 
(0.72 x 10 22 cirT 2 ). On the south-east of the remnant, region 12 
is faint, has a weak absorption (0.50x 10 22 cirT 2 ) and has a harder 
spectrum (2.17). 



bins of some X-ray profiles, there is not enough coverage to es- 
timate the flux. Those bins are removed from the profile. As in 
the spectral study, the two brightest point-like sources (see Sect. 
I3.21 i were removed from the X-ray images. The resulting radial 
profiles are shown in Fig.fTUland discussed in Sect. 14.3.21 

4. Results 

4.1. Radio flux 

An accurate estimate of the integrated radio flux density is crit- 
ical to investigate the mechanisms responsible for the high- 
energy emission. In this work we re-analyze the available radio 
data at 1.4 GHz following two different paths to estimate the 
radio flux density of RX J17 13.7-3946. 

Since the largest well imaged structure for the ATCA obser- 
vations is 25 arcmin in size, to recover information of structures 
larger than this size, the combination with single dish observa- 
tions is required. However, in this case, it was not possible to 
complete the procedure due to the fact that the overlap annulus 
of the interferometric and single antenna data seQin the uv space 
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1E-05 2E-05 3E-05 4E-05 5E-05 6E-05 7E-05 8E-05 



Fig. 4. Same X-ray mosaic as Fig.Q]smoothed to match the Point 
Spread Function of the HESS telescopes. The scale is linear 
and units are in ph/cm 2 /s/arcmin 2 . Overlaid are the y-ray con- 
tour excess (spaced at 30, 60 and 90 counts) from Fig. 7 of 
lAharonian et al.l (I2006I) . The sectors used for the radial profile 
comparison of Fig.[l0]are also drawn. 

is too small to produce a reliable image. We therefore used the 
single dish data at 1 .4 GHz to estimate the integrated flux den- 
sity within the area covered by the X-ray emission associated 
with the SNR (outer contour depicted in Fig. |5). This estimate 
provides an upper limit for the flux density of about 26 Jy with an 
uncertainty of the order of 10%. In addition we integrated over 
an interferometric image based on ATCA data within the same 
outer contour. In this case, since the largest well imaged struc- 
ture at this frequency is about 25 arcmin, we corrected "by hand" 
adding the minimum flux density required to fill in the few nega- 
tives remaining in the image. The surface brightness correspond- 
ing to this addition was below 7.4 X 10~ 24 W m~ 2 Hz-1 sr _1 . 
The off-source rms-noise is about 0.7 mJy/beam and the con- 
stant value added to fill in the few negatives is 1.2 mJy/beam. 
Such correction amounts less than 1 % of the total estimated flux 
density. 

Whether the radio emission from Arc 2 (see Fig. [5} is 
thermal or non-thermal in nature is an important issue as it 
is one the brightest features within the X-ray contours. In or- 
der to investigate the connection of Arc 2, with the SNR, we 
used mid-infrared observations. The right-hand panel of Fig. 
shows Spitzer 8 pm mid-infrared emission in the direction of 
RX J1713. 7-3946. Particularly, it can be seen intense infrared 
emission at the location of the HII region G347. 61 1+0.204. The 
infrared emission also evidences good morphological correla- 
tion with the radio Arc 2 suggesting a thermal origin for this 
feature. To investigate the nature of Arc 2 we applied the color- 
color criteria proposed by iReach et"ai] (120061) based on Spitzer 
data obtained at 3.6, 4.5, 5.8, and 8 pm, finding that Arc 2 has 
color characteristics compatible with polycyclic aromatic hydro- 
carbons (PAHs) origin (see Fig. Besides we applied HI ab- 



1 The only public single dish data available at this frequency come 
from the Survey of the South Celestial Hemisphere carried out with 



the 30 m dish o f the Argentine Institute of Radio Astronomy (IAR 
iTestori et alj200lb (HPBW-34 arcmin at 1 .4 GHz) 
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Fig.5. Left : Radio image of the region of RX J1713. 7-3946 at 1.4 GHz (from Lazen dic et al.l l2004). The scale is square root and 
units are in Jy beam 1 . Right : Spitzer image at 8 pm from the GLIMPSE survey. The X-ray contours are represented in black in 
both images and the radio contours in white. The regions used for the color-color plot of Fig.|6]are labeled here. The scale is square 
root and units are in MJy sr _1 . 



sorption techniques based on data extracted from the Southern 
Galactic Plane Survey (SGPS; iMcClure-Griffiths et alJ [2005b 
together with a flat rotation model for our Galaxy (assuming as 
solar parameters R e = 7.6 ± 0.3 kpc and = 214 ± 7 kms~'). 
From this study we conclude that the most distant HI absorp- 
tion feature for the Arc 2 is at -120 km/s, which corresponds to 
a near distance of ~ 6.7 kpc, placing in principle this thermal 
arc beyond RX J1713. 7-3946. With this information, to estimate 
the associated radio flux density we subtracted the contribution 
from all overlapping radio point sources (likely to be extragalac- 
tic and/or compact HII regions) and from Arc 2 (about 1.5 Jy) 
because it is likely to be unrelated with the remnant, estimat- 
ing a lower limit for the total flux density of ~ 22 Jy. Therefore, 
we conclude that the associated radio flux density at 1 .4 GHz is 
between 22 and 26 Jy. 



We note that th e 1-10 keV non absorbed flux obtained by 
iTanaka et all (l2008h with Suzaku is 47% higher (7.65 XlO" 10 erg 
cm" 2 s _1 ) than what we have found. However the absorption and 
index of their spectrum (0.79 x 10 22 cirT 2 and 2.39 respectively) 
are in agreement with our parameters. It is important to note 
that whereas our spectrum is extracted directly on the whole 
remnant, the Suzaku spectrum is the sum of spectra from 10 
particular regions scaled up to the whole remnant assuming 
the surface bright ness from the ASCA image (Sect. 3.3 of 
ITanaka et alJ [2008). In our spectrum, the point sources were 
removed but their contribution to the total flux is weak (less than 
1% for the Central Compact Object). We have cross checked 
the value of our absorbed global flux derived from the spectrum 
(in the 0.5-4.5 keV energy band) to the one derived from our 
mosaiced image. Both fluxes agree within 5%. 



4.2. Global X-ray flux 

We extracted a MOS spectrum from the whole remnant with- 
out degrading the spatial resolution of the data. The best-fit 
parameters (with background subtraction as described in Sect. 
[131 are N H = 0.66 ± 0.01 x 10 22 cnr 2 , index=2.37 + 0.01, an 
absorbed 1-10 keV flux of 3.95 ±0.03 x 10 -10 erg cnrV 1 and 
a non absorbed flux in the same band of 5.23 +0.04 x 10~ 10 
erg cm _2 s _1 . Our main source of uncertainty is the absolute 
calibration of XMM-Newton which is known with a precision 
of 10% rather than the very small statistical error (less than 
1%). Another flux estimate is obtained with a srcut synchrotron 
model using the radio flux upper limit derived in Sect. 14.11 The 
non absorbed flux derived is 5.49 +0.05 x 10~ 10 erg cm _2 s _I , 
very similar to the one obtained with a power law model. In the 
srcut model, the radio spectral index was fixed to 0.6 and a break 
frequency of 4.2 xlO 17 Hz was derived. This value, averaged 
across the remnant, is similar to the ma ximum break frequenc y 
reached in the bright limbs of SN 1006 dRothenflug et al.l l2004). 



4.3. X- and y-ray comparison 
4.3.1. Spectral results 

The best-fit parameters of the X-ray spectral modeling of the 
14 regions are given in Table [2] The large variation of photon 
index (1.9 < T < 2.6) seen in X-rays when using small extrac- 
tion regions in ICassam-Chenai' et alJ d2004l) have largely been 
reduced here with larger extraction regions and a degraded spa- 
tial resolution for the X-ray data (2.2 < T < 2.4, see Fig. [3). The 
comparison of the X- and y-ray photon index (Fig. [8} shows no 
significant correlation . The distribution of the photon index in 
X- and y-rays has a mean value of 2.32 and 2.09 respectively and 
a standard deviation of 0.075 and 0.073. Whereas the dispersion 
of the photon index at both energies is the same, the X-ray index 
is slightly higher than the y-ray one. However, there is a system- 
atic error on the y-ray photon index of 0.1 (AH06) that is to be 
added to the Fig[8] Whereas the variations of the photon index 
are small, there are significant variations in the flux from the 14 
large regions. Fig. [9] shows a good correlation between the X- 
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Fig. 7. Correlation plot between the X-ray photon index and the 
X-ray integrated flux (in the 1-10 keV band). The parameters 
shown here are the average of MOS 1 and MOS2 best-fit param- 
eters. The label over each point corresponds to HESS regions as 
defined in Fig. [2] 



Fig. 6. Infrared color-color diagra m for the regions defined in 
Fig.\5\(Right) overlaid on Fig. 2 of lReach et al.ld20 0r|). Two dis- 
tincts groups can be seen. The emission from the regions 8 to 12 
is compatible with PAH origin including our region of interest 
Arc 2 (labeled as region 9 here). Therefore the radio emission 
from Arc 2 is likely to be thermal and not related with the rem- 
nant. 

Table 2. Best-fit X-ray parameters obtained with an absorbed 
power law model for the HESS regions defined in Fig. [2] The 
X-ray flux is integrated in the 1-10 keV band. 



Region 


N H 


Index 


Integrated flux^ 




10 22 crrT 2 




(lO^crrrV 1 ~ 


1 


0.79 


2.43 (2.42-2.44) 


0.91 (0.90-0.92) 


2 


0.62 


2.35 (2.34-2.35) 


1.58 (1.57-1.59) 


3 


0.72 


2.35 (2.35-2.36) 


1.74(1.73-1.75) 


4 


0.76 


2.37 (2.37-2.39) 


0.76 (0.75-0.77) 


5 


0.59 


2.32 (2.31-2.32) 


1.42(1.41-1.43) 


6 


0.56 


2.44 (2.43-2.44) 


1.44(1.43-1.45) 


7 


0.67 


2.41 (2.40-2.42) 


0.89 (0.87-0.90) 


8 


0.64 


2.23 (2.22-2.24) 


0.58 (0.57-0.58) 


9 


0.48 


2.22 (2.22-2.23) 


0.90 (0.89-0.91) 


10 


0.46 


2.32 (2.31-2.33) 


1.15 (1.14-1.16) 


11 


0.63 


2.35 (2.34-2.36) 


0.72 (0.72-0.74) 


12 


0.50 


2.17 (2.15-2.18) 


0.51 (0.50-0.52) 


13 


0.51 


2.31 (2.30-2.32) 


0.69 (0.68-0.70) 


14 


0.57 


2.28 (2.26-2.29) 


0.70 (0.68-0.70) 



" Corrected for absorption 



ray integrated flux (1-10 keV band) and the y-ray integrated flux 
(1-10 TeV band). We note that the bright regions are brighter in 
X- rays than in y - rays. Such a behaviour has also been observed 
by iTanaka et al.l (120081) with Suzaku. However whereas they in- 
terpreted this as a linear correlation with some X-ray deviation, 
we interpret this as a non-linear correlation. We measured this 
non-linear correlation placing the X-ray flux in the X axis (as 
the error bars on the X-ray flux are smaller than those in y-rays). 
This measured slope is then inverted to obtain d log F*/d log Fy . 
The best-fit power law function F* = 0.81 + 0.09 (Fl) 2Al±0 - 55 is 




1.8 1.9 2.0 2.1 2.2 2.3 2.4 2.5 
Gamma-ray photon index 

Fig. 8. X-ray photon index against y-ray photon index. A sys- 
tematic error of 0.1 is to be added to the y-ray photon index. The 
dashed line is the bisector. 



statistically preferred (x 2 /d.o.f. = 9.28/12) over the best-fit lin- 
ear function F* = 1.01 ± 0.04 (x 1 1 d.o.f. = 61.06/13). We 
note that the fit with a linear function Ff = aF y v - b (bestfit pa- 
rameters : a = 2.26 and b — 1 .43) gives similar results in term of 
X 1 (x 2 /d.o.f . = 9.38/12) than the non-linear fit but is not physi- 
cally understandable since both images are already background 
subtracted. 



4.3.2. Morphological results 

The comparison of the radial profiles allows us to investigate the 
extent of each emission as well as to localize their respective 
peaks. Such a study has already been carried out by AH06 by 
comparing the ASCA and HESS data. However as the coverage 
with ASCA did not always reach the boundaries of the SNR it 
was not possible to compare the extension of the remnant in both 
wavelengths. 
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Fig. 9. Log-log correlation plot between the HESS y-ray inte- 
grated flux (1-10 TeV band) and the XMM-Newton X-ray in- 
tegrated flux (1-10 keV band, using MOS1&2 data). We can 
see that the correlation is best-fitted with a non-linear function 
= 0.81 ± 0.09 (Fl) 2A1±0 - 55 (solid line). The best-fit linear 
function = 1.01 ± 0.04 F y v is also represented (dashed line). 



The global agreement in the 8 sectors is good and particu- 
larly striking in sector 7 (the brightest spot of the remnant in 
both wavelength) as can be seen in Fig. [10] In regions 5 and 
6 we see that the peak of the emission, along the radial direc- 
tion, is located closer to the center in X-rays than in y-rays. The 
global shape of the radial profile is similar in both wavelength 
but shifted by ~0.1° to the center in X-rays. 

It is important to note that we did subtract an astrophysical 
background to the X-ray radial profiles. To estimate the back- 
ground, we used the same region as in the mosaic building sec- 
tion (Sect. [3J)- 



5. Discussion 

5.1. Global synchrotron spectrum 

The radio flux density at 1 .4 GHz that we derive in our study (22 
Jy < S < 26 Jy) for the whole remnant is significantly higher 
than the flux density previously published . However it is im- 
portant to note that the flux calculated by Elli son et alJ (1200 ll) 
(S=4+l Jy) was estimated only locall y for the two bright arcs 
in the northwestern region. Later on, lLazendic et alJ d2004l) de- 
rived for these features a flux density of 6.7 ± 2.0 Jy based on 
the improved radio image at 1 .4 GHz. In X-rays the contribution 
of the northwestern region accounts for ~ 20 % of the total flux 
in the 0.5-4.5 keV energy band. If the ratio of synchrotron emis- 
sion in X-rays to that in radio is more or less constant throughout 
the remnant, the radio flux expected for the whole remnant is of 
the order of 26 J>0 which is compatible with our estimation. We 
note that the flux that we derive is about 2 times higher than the 
flux of 13.4 Jy used in AH06. This flux density was obtained 
assuming that the flux density for the en tire remnant wa s abou t 
twice the flux density value obtained by lLazendic et al. (2004) 



2 Without taking into account the thermal emission from Arc 2 
that accounts for ab out 1.5 Jy out of the 6.7 Jy estimation from 
lLazendic etal] d2004f) . 



for the northwestern region. Our higher new estimate tend to re- 
duce the difference between the measurements and predictions 
of the h adronic models usua lly requiring high magnetic fields > 
100 nG dTanaka et al.ll2008b . 

In X-rays, the 1-10 keV non absorbed flux that we derive with 
XM M-Newto n is sig nificantly lower than the value derived by 
iTanaka et al] (l2008h with Suzaku (5.23 and 7.65 xl0~ 10 erg 
cm" 2 s _1 respectively, see Sect. 14.21 . This new X-ray flux to- 
gether with the higher radio flux estimate decrease the global 
X-ray to radio ratio by about a factor of 3. Such a change is 
clearly important in the shape of the synchrotron spectral energy 
distribution and could impact the results from multi-wavelength 
models of the remnant. 



5.2. X- and y-ray comparison 
5.2.1. Spectral index 

When using large extraction regions and degrading the spatial 
resolution of the X-ray data, we found no significant spatial vari- 
ation of the photon index (see Sect. I4.3.1I I. In other words, the 
small-scale variations reported by ICassam-Chenai et alJ (2004) 
are washed out at the resolution of the HESS telescope. If there 
exists small-scale spectral variations in y-rays, one expect them 
to be washed out as well with the current HESS resolution. 
However the mean X-ray photon index (2.32) is slightly steeper 
than the y-ray one (2.09). In the case of IC off the CMB photons 
for the y-ray emission, the mean electron energy required to have 
a photon at 1 TeV is Ef-= 16 TeV (for IC off IR or optical pho- 
tons it is of course lower). The electrons emitting synchrotron at 
1 keV with a magnetic field of 70 fiG have an energy of E S J DC = 
27 TeV. At those energies the electrons are close to the cutoff and 
we do expect a higher index when looking at electrons of higher 
energies. We can reasonably set an upper limit to the magnetic 
field as El ync cannot be lower than E l J~ in order to reproduce the 
steeper index of the synchrotron spectra. This sets an upper limit 
on the downstream magnetic field of 200 //G, consistent with the 
estimate coming from the width of the X-ray fi laments (about 70 
/jG from lBallefcOOrjllBerezhko & Volkll2006h . 



5.2.2. Correlation between the X-ray and y-ray fluxes. 

We have found in Sect. 14. 3.T1 that the X-ray flux is indeed corre- 
lated with the y-ray flux, but the correlation appears non-linear, 
more like dlogF^/dlogF? = 2.41 ± 0.55 (Fig.©. We discuss 
here why it could be so. We will assume that most of the range in 
flux that we see is due to density variations around the remnant. 
It is certainly not entirely true (geometric effects such as limb 
brightening must play a role at some point) but it is probably 
indeed the major contributor in view of the very structured in- 
terstellar gas and X-ray image. We will assume that in the range 
which emits the y and X-ray (i.e. one to a few hundred TeV) 
the particle distribution can be represented as a cut-off power 
law dN/dE = KE S C(-E /E c ) in which K is proportional to the 
ambient density and F c is limited by synchrotron losses for elec- 
trons (£ c , e oc y sh / y[Bi) and by age for protons (F c , p oc BdV^ h ) 
where Z?d is the downstream magnetic field. The shape of the 
cutoff C(x) can be any function decreasing from 1 to 0. We will 
further assume that magnetic field may increase with density 
(fid oc yfi) and that the shock velocity adjusts as V s h oc «~ - 5 
(same pressure). 
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Fig. 10. Radial profiles in X-rays in two energy bands and in y-rays for the 8 sectors defined in Fig. [4] The general agreement is 
good and particularly striking in sector 7 (the brightest spot of the remnant in both wavelengths). However there are also interesting 
differences in sectors 5 and 6 where the bulk of the X-ray emission seems to come more from the inside of the SNR than in y-rays. 



By calculating dlogF v /dlog« for the synchrotron and 
Inverse Compton case (see Appendix lAl. we can then predict 
in the flux-flux correlation comparable to Fig. [9]: 

dlogfT C _ ¥(l+P)-ax 
dlogFje 1 - (1 +p)(a y - £=1) 

where ax and a y are respectively the X and y-ray spectral slope 
in energy not photons (a, = -d log F' v /d log v). 

For a standard value of s — 2 and the observed values of ax 
= 1.32 and a y = 1.09 this gives (0.18 + 1.5 B) / (0.41 - 0.59 B). 
If B& is insensitive to density (B = 0) the predicted correlation 
is opposite of what is seen: the range in X-ray flux would be 
smaller than in y-rays, because the negative feedback on v c via 
V s h plays more strongly in X-rays which are further in the cutoff 
part of the spectrum. But a modest dependence of Z?d on density 
like B = 0.1 is enough to invert the trend because in the loss- 
dominated regime the y-ray cutoff frequency decreases with Z?d 
while the X-ray one is independent of B&. In other words, the 
slope of the \og(Fx) vs \og(F y ) correlation is very sensitive to B. 
To get the observed value of 2.41 requires B = 0.28. 

In view of the oversimplified character of that approach we 
do not claim that this is a measurement of dlogBd/dlogn but 
we think it shows that such a steep correlation is reasonable in 
a leptonic model. The specific model in which we have tried to 
push the calculation further (Appendix IB1 does not give a con- 
sistent answer, but it is far from unique. One way to improve on 
the measurements would be to use extraction areas in which the 
filling factor of the SNR is the same, like the angular sectors on 
Fig HI This would leave in the flux variations only what is due 
to varying external conditions. It requires reanalyzing the HESS 
data, so it is left for future work. 

If we now turn to the hadronic hypothesis, the same line of 
reasoning (see AppendixlAli then leads to : 

d l 0g /rhadr 2-(\-B)(a y + \-s) 

For a standard value of s — 2 and the observed values of ax 
= 1.32 and a y = 1.09 this gives (0.18 + 1.5 B) / (1.91 + 0.09 B). 



It is clear that whatever/? < 1 (it is hard to imagine how B^ could 
increase faster than density) this quantity is always < 1 . In other 
words, no magnetic field can make up for the natural n 2 charac- 
ter of the hadronic mechanism which predicts a fast increase of 
^■hadi w j t j 1 density, go a t i eas t m th a t (over)simplified framework 
the correlation we observe is not in favor of a hadronic model. 

There remains the possibility (Malkovetal. 2005) that the 
density increases very fast outwards (SNR hitting a shell) to 
the point where most of the y-ray emission arises outside the 
remnant (in the precursor). In that case the width of the pre- 
cursor increases as E so that the spectral shape of the y- 
ray emission may be estimate d by multiplying Fy aAl ~ by E 
( Zirakashvili & Aharonian 2007). This in turn changes (O into 

djogFT _ ^(l+B)-a x 
dlogFj! adr 2-(\-B){a y + 2- s) 

so that for s — 2 and the observed spectral indices one expects 
(0.18 + 1.5 B) I (0.91 + 1.09 B) for the slope of the log(F z ) vs 
\og{F y ) correlation. This is still always < 1 . 

5.3. Spatial comparison 

The main difficulty of a leptonic model to account for the obser- 
vations in RX J1713. 7-3946 is that it requires a small magnetic 
field (on the order of 10 fiG ; AH06) to explain the rather large 
y to X-ray ratio, if the emitting volume is the same. This is in- 
consistent with the magnetic field derived from the width of the 
X-ray filament s (70 a <G or so). A possible reason, suggested by 
lLazendic e t al. (2004), is that the magnetic turbulence decays be- 
hind the shock faster than the electrons lose energy. This leaves 
a larger volume (downstream) to IC than synchrotron, and does 
not require that large a magnetic field to begin with. A definite 
prediction is then that the y-ray emission should peak inside the 
X-rays. 

In the comparison of the X- and y-ray radial profiles (14.3.2b 
we did see a radial shift, particularly visible in region 6 (West), 
between the X- and y-ray emission. But the shift is in the op- 
posite direction, i.e. the X-ray emission peaks at smaller radius 
than the y-ray emission. The value of this shift for region 6 is 
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~ 0.06°=1.2 pc with a remnant at 1 kpc. Another effect that 
we have not discussed is that because of the magnetic jump at 
the shock (typically a factor 3 if the magnetic field is mostly 
turbulent and isotropic), the synchrotron emission is strongly 
suppressed ahead of the shock whereas the IC emission will 
decrease more smoothly over one diffusion length. In a Bohm 
regime the diffusion coefficient is D% = r^c/3 where is the 
Larmor radius and the diffusion length Zdiff = D^/V^ where V s h 
is the shock wave velocity. For an electron energy of 16 TeV, a 
shock speed of 4000 km/s and an upstream magnetic field of 10 
/iG we have Zdiff = 0.1 pc. This is small in comparison with the 
shift of 1.2 pc. 

However there exists another purely geometric possibility to 
explain why the X-rays peak inside the y-rays. Actually from 
FigUHlthe effect is significant only in regions 6 (West) and 8 
(North). In both regions it looks from the image (FigQ]) that the 
remnant extends beyond the main X-ray peak. This is typical 
of a line of sight superposition of a region with larger density 
(brighter, slower shock) and a region of lower density (fainter, 
faster s hock). This idea is supp orted by the fact that observations 
in CO dMoriguchi et al.ll2005l) show that clouds are present in 
those particular regions. Now from Fig(7]we know that the X- 
ray brightness increases much faster than the y-ray brightness. 
The X-ray radial profile is then dominated by the inner bright 
edge, whereas the outer plateau may contribute significantly to 
the y-ray emission and shift the integral along the line of sight 
to a larger radius. Because of that possibility we cannot say that 
the radial profile is incompatible with a leptonic model. 

In a hadronic model the simplest way to understand that 
the y-ray emission is further out is if a sizable fraction of the 
flux come s from outside the sh ock. This is possible and actually 
expected (Mal kov et al.l l2005h if the ambient density increases 
ahead of the shock (SNR entering a molecular cloud). There are 
then vast amounts of upstream gas to be used as targets for ac- 
celerated protons. A definite prediction of that model is that the 
remnant should look larger as energy increases because it will be 
governed by the diffusion length. The energy range in the HESS 
data may be enough to test that, and certainly will after HESS 2. 

6. Summary 

New XMM-Newton observations of RX J1713. 7-3946 have al- 
lowed us to complete the coverage of the remnant. We have then 
carried out a comparison of the remnant in X- and y-rays. The 
radio emission has also been studied. This leads us to the follow- 
ing conclusions : 

1. The small scale variations of the X-ray photon index are 
largely smeared out at the resolution of the HESS telescope. 
This is consistent with the non detection of photon index 
variation in the y-rays. 

2. The mean X-ray photon index (2.32) is slightly steeper than 
the y-ray one (2.09). This is expected in a leptonic scenario 
as the electrons emitting y-ray Inverse Compton have a lower 
energy than those emitting X-ray synchrotron. The photon 
index offset is thus the result of the spectrum cutoff of the 
electrons. 

3. The comparison of the X-ray vs y-ray integrated flux sug- 
gests a non-linear correlation with dlogF*/dlogFy = 2.41 
+ 0.55. If the range of flux that we see is due to the variation 
of the density around the remnant then a leptonic model can 
more easily reproduce the observed X/y-ray correlation. 

4. The comparison of the radial profile in X- and y-rays indi- 
cates that for some regions, the X-ray emission comes more 



from the inside of the remnant than in y-rays. This can be 
explained in a hadronic model if a fraction of the y-ray flux 
comes from the outside of the remnant (interaction of the 
remnant with a cloud). However this radial shift could also 
be due to superposition effects. Therefore we can not rule out 
the leptonic scenario. 
5. Concerning the radio counterpart of the remnant, we have 
shown that one of the brightest arcs seen at 1 .4 GHz (Arc 2) 
is of thermal origin and likely not associated with the rem- 
nant. Taking this into account, we have derived a lower limit 
on the integrated flux density of 22 Jy and an upper limit of 
26 Jy. 
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Appendix A: Correlation between the X-ray and 
y-ray fluxes 

The synchrotron emission can be written 

Fr*K e B^v-^s[-^ (A.1) 

where S is a function characteristic of synchrotron emission and 
the shape of the electron cutoff and v s J nc oc BE\ e cc V^ h . Then 

diog/y 3 ™ _ s+l _ 5^ _v_ dlogVe ync 
dlogn 2 P S v7 nC dlogn 1 ' ; 
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Noting that 

diog£T;_ 

dlog v 



s- 1 S' 
~~2~ + J v; 



sync 



= -ax 



(A3) 



is the X-ray spectral slope (in energy, not photons) and that 
d log v s J nc I d log n = - 1 , we can write 

dlog/v 5 " 10 s + 1 s-1 s + 1 

■ f = 1 + -^-P ~ »x + = -=-(1 + B) - a x (AA) 
dlogn 2 2 2 

As long as the Klein-Nishina reduction of the cross-section 
is not reached (this is still true for 10 TeV electrons on the CMB) 
the Inverse Compton emission may be written in the same way 
as ( IA. 1 b but without any explicit B term and with v' c oc E 2 e oc 
V 2 ./B. The same line of reasoning then leads to 



dlogn 



ic 



= l-(\+B)\a 7 - 



(A.5) 



in which a y is the y-ray spectral slope. 

Putting together (lA.4b and (IA.5t we then predict in the flux- 
flux correlation comparable to Fig(9] 



sync 



d log F' Y 



^-(\+B)-a x 



dlo g n c 1 - (1 +B)(a y - *=1) 



(A.6) 



If we now turn to the hadronic hypothesis, the n° decay emis- 
sion may be written 



lyhadr 



(A.7) 



where H is a function characteristic of n° emission and the shape 
of the proton cutoff and 
reasoning then leads to 



of the proton cutoff and v£ adr oc F cp oc BV 2 h . The same line of 



dlogFjT 
d log n 



= 2-(l -B)(a 7 + l-i) 



(A.8) 



in which we now assume that hadronic emission dominates the 
y-ray spectral slope. Putting together (|A.41 > and (|A. 8I > we then 
predict for the flux-flux correlation in the hadronic model 



dlogF; 
dlogF* 11 



sync 



!f(\+B)-ax 
2-(\-P)(a y + \-s) 



(A.9) 



There remains the possibility (Mal kov et al.ll2005h that the 
density increases very fast outwards (SNR hitting a shell) to 
the point where most of the y-ray emission arises outside 
the remnant (in the precursor). In that case the width of the 
precursor increases as E so that the spectral shape of the 
y-ray emission may be estima ted by multiplying IA.7I by E 
( Zirakashvili & Aharonian 2007). 

This in turn changes ( |A.8t into 



(A. 10) 



= 2-(l -B)(a y +2-s) 



dlog 
dlogn 

and changes ( |A.9t into 

dlogFr c = f(l+B)-a x 
dlogF^ dl " ~ 



2-(l -B)(a y + 2-s) 



(A.ll) 



Appendix B: Detailed calculation for electrons 
dominated by radiative losses 

If the cutoff in the electron spectrum is defined by synchrotron 
cooling as we assume in Sect |5.2.2"l and ApplAl then the cooling 
continues downstream and the electron distribution integrated 
over space gets steeper (the power law index increases by 1) 
down to a break energy E\, oc B' 2 ^ 1 in which fo is the SNR age 
dZirakashvili & Aharonianll2007l) . Since we are interested only 
in what happens above E\, this can be accounted for in ( IA. 1 b by 
changing s to the steeper value and considering that the normal- 
ization K e follows £b (the total number of electrons is always 
dominated by those below Fb) or B' 2 (to is the same in all parts 
of the SNR). This amounts to adding -2B to Eqs dA.4| ) and ( |A.5t , 
resulting in 



d log F 



sync 



S+ 1 



dlogn 



s-3 
a x + —z-B 



ic 



d log F; 
d log n 



5+1 



B ar r + 



For the most interesting case j = 3we then have 



sync 



dlogF, 

ic 



2 - a x 



2 - a y - B( a y + 1) 



(B.l) 
(B.2) 

(B.3) 



d log F 

Fora x = 1.32 and a y = 1.09 we get 0.68/(0.91 - 2.09/3). This is 
very different from what is obtained with s - 2 and no break in 
the electron spectrum, but reaches the observed value 2.41 for a 
very similar value of B — 0.30. 

The pion decay emission is insensitive to syn- 
chrotron cooling so that "radiative" model would predict 
dlogF' ync /dlogF^ dr = 0.68/(1.91 + 0.09/?). This is again 
always less than 1 . 

The same framework also naturally predicts variations of the 
spectral index going with the flux variations. Computing that re- 
quires a specific representation of th e spectral shape. Suitable 
appro ximate formulae are given by [Zirakashvili & Ahar onianl 
d2007l) in their Eqs (35) and (46). They apply when there is no 
jump in magnetic field at the shock. This is not the most likely 
situation in our opinion, but it is still interesting to carry out the 
exercise to the end. The observed average spectral slope in X- 
rays ax = 1.32 corre sponds to x = v/v?™ = 4.5, in keeping with 
what was derived by lTanaka etaP d2008l) (v' ync = 0.67 keV). In 
the same way, the observed spectral slope in y-rays a y = 1.09 
corresponds to y - v/v^. c = 4. 1. From Eqs (34) and (45) of 
IZirakashvili & Aharonianl d2007l) we note that v^/vj nc = (1.2 
TeV) / (2.2 keV) (B/100/iG)" 1 . For B 70juG and vj nc = 0.67 
keV this predicts v' c - 0.5 TeV which is reasonable. 

Coming back to the spectral index variations, since 
d log Vj ync /d log n = — 1 when it is limited by cooling, we have 



dffx 
d log n 



dcx d log v c y da 

= x 

d log x d log n dx 



(B.4) 



0.25 jc ' 5 - 0.3795 x 6 (l + 0.46 x 6 )' 



(B.5) 



dcrx 
d log n 

At x = 3.5, this gives dcx/dlogn = 0.32. It then predicts 
d log Fy yac /dax - 0.74 / 0.32 = 2.3. The peak to peak dispersion 
that is observed on the X-ray flux A log F v yac 1 .2 should then 
be associated with a peak to peak dispersion on the slope Aax = 
0.53. This is larger than what is observed on Fig.^Acry = 0.27). 
This is qualitatively expected because flux variations are also af- 
fected by the geometry (the regions do not cover the same frac- 
tion of the SNR). However the small amplitude of the index vari- 
ations implies that the fraction of the flux variations that is due 
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to geometry A log Fy eom is large. It should add quadratically with 
the part due to density variations (2.3 x Aa x ) so A log Fy eom ~ 
1 . 17. Of course geometrical effects should affect the y-ray emis- 
sion in the same way but the observed dispersion on y-ray flux 
is only A log F v ync 0.6. This means we are at a dead-end. That 
specific model cannot explain at the same time the relatively uni- 
form X-ray spectra and the larger contrast in the X-ray flux than 
the y-ray flux. 
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Fig. 1. EPIC MOS plus PN image in the 0.5-4.5 keV band. The units are ph/cm 2 /s/arcmin 2 and the scale is square root. The image 
was adaptively smoothed to a signal-to-noise ratio of 10. 



